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Abstract: We present a photometric study of the neglected open cluster Berkeley 53. We derived its fim- 
damental parameters, such as the age, the interstellar reddening, and the distance from the Sun, based on BV 
photometry combined with near-infrared JHK^ data. The structure and the mass function of the cluster were 
also studied and the total number of members and the total mass were estimated. The cluster was found to be 
a rich and massive stellar system, located in the Perseus Arm of the Milky Way, 3.1 ± 0.1 kpc from the Sun. Its 
age exceeds 1 Gy but it seems to be very young in the context of its dynamical evolution. The analysis of the 
two-color diagrams and color-magnitude diagrams indicates that the cluster is significantly reddened. However, 
both methods resulted in different values of E{B — V), i.e. 1.21 ± 0.04 and 1.52 ± 0.01, respectively. This 
discrepancy suggests the presence of an abnormal interstellar extinction law toward the cluster. 
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1 Introduction 

The open cluster Berkeley 53 (C 2055-1-508) was discovered by Setteducati & Weaver (1960). Ruprecht (1966) 
classified it as a poor, concentrated open cluster of Trumpler type II3p. According to the New catalogue of 
optically visible open clusters and candidates by Dias et al. (2002), the cluster's apparent diameter is 12' and 
Trumpler type is III2m. No dedicated studies of this object have been performed to date. 

Berkeley 53 consists of stars fainter than V = 18 mag and is located in the vicinity (3.'1 from the cluster 
center) of the bright {V = 6.6 mag) foreground star HD 199578. This makes photometric observations of the 
cluster difficult. 

In this paper we present a photometric study of Berkeley 53 resulting in the determination of its basic 
parameters, such as the age, the interstellar reddening, and the distance from the Sun. The structure and the 
mass function of the cluster are also studied and the total number of members and the total mass are estimated. 



*This work is based on observations made with the 2-m telescope of the Rozhen National Astronomical Observatory, which is 
operated by the Institute of Astronomy, Bulgarian Academy of Sciences. 
tgm@astri.uni.torun.pl 
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2 Observation and data reduction 



Observations were performed with the 2/16 m Ritchey-Chretien telescope of the Rozhen National Astronomical 
Observatory (NAO, Bulgaria), operated by the Institute of Astronomy, Bulgarian Academy of Sciences. The 
instrument was used in a direct imaging mode and was equipped with a Princeton Instrmnents VersArray:1300B 
CCD camera mounted in the Ritchey-Chretien focus. The field of view was 5.'8 x 5.'6 with a scale of 0.26 arcsec 
per pixel. Observations were carried out on August 19, 2007. Two exposures in the B filter and three frames in 
the V band were acquired. The exposure time was 600 s and 300 s, respectively. Four exposures in the U filter 
with 900 s exposure time were also acquired but they occurred to be not deep enough to detect the cluster stars. 

The field around Berkeley 53 is presented in Fig. [T] where the fragment of the sky covered by our observations 
is also sketched. The telescope was not pointed on the cluster center due to the nearby bright star HD 199578. 
The B band exposures were affected by a narrow strip of its reflected light passing horizontally through the 
center of frames. 

CCD frames were processed using a standard procedure that included subtraction of bias frame, fiat-fielding 
with twilight flats, aperture and point spread function (PSF) photometry, transformation to the standard system, 
and astrometric calibration. The aperture photometry and astrometric calibration were performed with the 
software pipeline developed for the Semi- Automatic Variability Search sky survey (Niedzielski et al. 2003). 
Exposures in a given band were averaged. The PSF magnitudes were obtained with IRAF0 package DAOPHOT. 
A second-order variable PSF was used to compensate for variability of stellar profiles across frames. Aperture 
corrections were determined using aperture photometry of 13 isolated stars that were used while building the 
PSF profile. 

The raw instrumental magnitudes were corrected for the scattered light effect (Markov 2008) . The following 
correction formulae were obtained by us: 

A6= (0.122 ±0.006) (1) 

Au = (0.106 ±0.004)^2, (2) 

where A6, Aw are the additive magnitude corrections in the corresponding filters and d is the distance of a star 
from the CCD matrix center normalized by the detector half-size. 

The calibration coefficients that transform instrumental magnitudes into standard ones were determined 
using 16 stars from the Messier 92 field (Majewski et al. 1994, see their Table 1) and 40 stars from the field 
of NGC 7790 (Odewahn et al. 1992). The magnitude range was between 13.144 and 18.315 mag in V and the 
[B — V) color index coverage was in the range between —0.111 and 1.918 mag. The following equations were 
derived: 

b-B = (0.58 ± 0.03)A: - (0.17 ± 0.01)(S -V)- 23.06 , (3) 

v-V = (0.33 ± 0.02)A: - (0.16 ± 0.01)(S -V)- 23.67 , (4) 

where 5, v are the corrected instrumental magnitudes, B, V are the magnitudes in the standard system, and X 
is the airmass. 

The final list of stars contains equatorial coordinates, V magnitudes, and [B — V) color indices. It is available 
in electronic form at the survey web sitfi and the WEBD^Ifl database (Mermilliod 1996). 

^IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation, 
^http: / /www.astri.uni.torun.pl/"gm/OCS 
'^http: / /www.univie.ac.at/webda/ 
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Figure 1: A 9' x 9' field around Berkeley 53. The 5/8 x 5.'6 field of view of the telescope is marked with a 
square frame. The redetermined center of the cluster is marked with a cross. The core radius boundary is also 
sketched. The bright star at the left edge of the image is HD 199578. The image was taken from the Digitized 
Sky Surveys (DSS). North is up, East to the left. 

The optical data were complemented with near-infrared JHKg photometry extracted from the 2-Micron All 
Sky Survey (2MASS, Strutskie et al. 2006). The extraction radius was set to 30' around the cluster center (see 
Sect. 3.1). Optical photometry was combined with near-infrared data to perform a comprehensive photometric 
study of the cluster (e.g. Maciejewski 2008). 



3 Data analysis and results 
3.1 Cluster structure 

The cluster structure was investigated with the radial density profile (RDP). For this purpose only the 2MASS 
data set was used for constructing RDP due to unlimited field of view. We started with redetermining the 
cluster center. The algorithm was adopted from Maciejewski & Niedzielski (2007) and started from a tentative 
position determined by eye. Two perpendicular stripes were cut along declination and right ascension starting 
from the approximate cluster center and the histogram of star counts was built along each stripe. The bin with 
the maximum value in both coordinates was taken as a new cluster center what resulted in RA = 20'^55™56^ and 
DEC = +51°02.'8 for epoch J2000.0 {I ^ 90?29, b = 3?75). As one can see in Fig. [H these values do not point 
to the center of the field of view. Next, the profiles were constructed by counting stars inside concentric rings 
of width 1', centered at the redetermined cluster center. The density uncertainty in each ring was estimated 
assuming Poisson statistics. The RDP is plotted in Fig. [51 

To parametrize the density distribution, a two-parameter King (1966) density profile was fitted with the least- 
squares method in which the uncertainties were used as weights. We derived the core radius (the distance where 
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Figure 2: The radial density profile. The solid line denotes the fitted density distribution while the dashed one 
marks the level of the stellar background density. 
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Figure 3: (H — Kg) vs. (B — V) and (J — H) vs. (H — Kg) diagrams constructed for the stars located within the 
core radius from the cluster center. The stars that form the main sequence of the cluster arc marked with filled 
circles. The main sequence, fitted by shifting along the reddening vector, is sketched with a continuous line in 
each diagram. 

the stellar density drops to half its maximum value) Tc = 2.5 ± 0.1 arcmin, the central density /o = 18.0 ± 0.5 

stars/arcmin^, and the background density level /bg = 6.9 ± 0.1 stars/arcmin^. The fitted profile is sketched 
with a solid line in Fig. 2 while the dashed line marks the level of /bg. The chistcr limiting radius was 
roughly estimated by eye-inspection in the RDP. The stellar density excess is visible up to at least 11' - a value 
almost two times greater than the literature one. 

The RDP allowed us to estimate number of observed stars belonging to the cluster to be -^900. That suggests 
that Berkeley 53 is a very rich and massive stellar system. 

3.2 Two-color diagrams 

{H — Ks) vs. {B — V) and (J — H) vs. {H — Ks) two color diagrams (TCDs) were constructed to estimate the 
interstellar reddening toward the cluster. Only stars located within the core radius from the cluster center were 
considered to minimize the influence of the background star contamination. Among them only stars forming the 
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main sequence of the cluster were taken for the further analysis. They are marked with filled circles in Fig. [3] 
where the TCDs are plotted. 

The theoretical main sequence was extracted from the Padova isochrones for solar metallicity Z = 0.019 
(Giraldi at al. 2002). The reddenings E{B-V) = 1.18 and E{H-Ks) = 0.26 were obtained by shifting the main 
sequence (a continuous line in Fig. 3) along the reddening vectors whose normal slopes were calculated assuming 
the universal interstellar extinction law by Schlegel et al. (1998). To obtain an independent determination of 
the color excess in {B — V) color index, the value of E{H — Ks) was transformed into E{B — V) applying the 
relation ^^(^lyf = 0.209 taken from Schlegel et al. (1998). We derived E'{B - V) = 1.24. The resuhs seem to 
be consistent with a mean value of {E{B — V)) = 1.21 ± 0.04. 

3.3 Color-magnitude diagrams 

A preliminary J vs. (J — Kg) color-magnitude diagram (CMD) for overall cluster region (r < r^^) is presented 
in Fig. m We applied a decontamination procedure to remove background-star contamination. The details of 
this procedure can be found in Maciejewski & Niedzielski (2007). The CMDs were built for the cluster region 
and for an offset field. A concentric offset field of width 10' and starting at r = rum + 2' from the cluster center 
was used. Then the CMDs were divided into two-dimensional bins and the number of stars within each box was 
counted. The cleaned (decontaminated) cluster CMD was built by subtracting the number of stars in the offset 
box from the number of stars in the corresponding cluster box. The latter number was weighted by the cluster 
to offset field area ratio. Knowing the number of cluster stars occupying each box, the algorithm randomly chose 
the required number of stars located in the cluster area and with the adequate magnitude and color index. 

As one can see in Fig. 01 a rich main sequence is clearly visible, as well as a red giant clump. The estimated 
number of observed cluster stars is ~1000 after rejecting stars with outstanding magnitudes and colors. This 
value is comparable with the one obtained from the RDP analysis. We found ~200 evolved stars constituting 
the red clump. 

To determine fundamental astrophysical parameters of the cluster, we built 8 CMDs combining V and J 
magnitudes with (B — V), (V — J), {V — H), and {V — Ks) color indices. The broad color baselines were expected 
to minimize photometric errors and to determine precisely the reddening. Only stars located within the core 
radius around the redetermined cluster center were considered to minimize the influence of the background star 
contamination on further analysis. The diagrams are plotted in Fig. [5l 

Initially, we tried to derive the distance modulus and the age of the cluster via fitting a set of theoretical 
Padova isochrones (Giraldi at al. 2002). Since there is no information about the cluster metallicity, the solar 
metallicity oi Z — 0.019 was assumed. The fitting algorithm is based on the least-squares method and uses stellar 
magnitudes as weights. Stars with extremal color indices or magnitudes were rejected manually before running 
the isochrone-fitting procedure. The fitting and decontamination procedures were run independently for each 
diagram with a fixed value of reddening. The latter quantity was calculated for a given color using {E{B — V)) 
from the two-color diagram analysis and assuming the universal interstellar extinction law by Schlegel et al. 
(1998). However, we were unable to find a convergent and satisfactory final solution and the reddening was 
found to be underestimated. 

To explain this discrepancy, we ran the isochrone-fitting algorithm with reddening as a third free parameter. 
The results of individual fits are collected in Table 1 and the respective isochrones are drawn with solid lines in 
Fig. 5. The dereddened distance modulus (to — M)o was calculated assuming the universal interstellar extinction 
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Figure 5: The color-magnitude diagrams constructed for the central (r < r^) part of Berkeley 53. The best-fit 
isochrone is drawn with a solid line and the zero age main sequence is sketched with a dashed one. 
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Table 1: The results of a solar-metallicity isochrone fit for individual CMDs. E denotes the fitted value of the 
color excess for a given color, (m — M) - the fitted value of the apparent distance modulus, E{B — V) - the 
calculated E{B — V) color excess, and (m — M)o - the calculated dereddened distance modulus. 
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law by Schlcgcl ct al. (1998) with a total-to-selcctivc absorption ratio of i? = 3.315. The value of E{B — V) was 
found to bo about 0.3 mag greater than the two-color diagram analysis suggested. One can note that E{B — V) 
tends to be slightly smaller for colors including more reddish bands. We obtained e(b-v) ~ 2-414 ± 0.025, 
E^Zy] = 2.707±0.017, and ^^Jb-v) = 2.893±0.023. Comparing to the values given by Schlegel et al. (1998), 
i.e. 2.413, 2.739, and 2.948, respectively, a discrepancy can be noted in bands H and and increasing with 
wavelength. This effect explains the lower value of the interstellar reddening obtained from the TCDs analysis. 
Adopting redetermined reddening vectors, wc obtained consistent results of TCDs and CMDs analysis. In the 
absence of photometric observations in the U band, it is impossible to answer the question about the nature 
of the observed effect. It may be caused by an abnormal interstellar extinction law toward Berkeley 53 or by 
systematics in 2MASS photometry. 

Finally, the following mean results were obtained from the CDMs analysis: log(age) = 9.09±0.02, E{B — V) = 
1.52 ± 0.01, (m — M)o = 12.45 ± 0.05, and a distance of 3.1 ± 0.1 kpc. The linear diameter was found to be 
19.8 ± 0.5 pc. The main possible cause of systematic error in the values given above is the unknown metallicity 
of the cluster. To estimate its influence on the results, we repeated the isochrone-fitting procedure for the super- 
and sub-solar metallicities oi Z = 0.030, Z = 0.008, and Z = 0.004. We noticed the increase of the reddening 
with the decrease of the metallicity {AE(B — V)^ 0.07 for each metallicity step, including Z = 0.019) while the 
age and the apparent distance modulus remained stable. It is also worth noting that the discrepancy between 
the TCDs and CMDs values of the reddening cannot be justified by non-solar metallicity of the cluster. 



3.4 Total mass and number of members 

Studies of the mass function (MF) were carried out to estimate the total mass of the cluster and the number of 
its members. The analysis is based on the 2MASS photometry due to the wide field of view available and uses 
the algorithm adopted from Maciejewski & Niedzielski (2007). The first step was to build a luminosity function 
(LF) for overall (r < rum) cluster region. The bright end of the LF was determined by the main-sequence 
turn-off point while the faint end was set for J = 15.8 mag - the value of the 99.9% Point Source Catalogue 
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Figure 6: The mass function derived for the main-sequence stars of Berkeley 53. 



Completeness Limii|f|. We used bins as small as 0.1 mag due to cluster richness. Another LF was built for a 
concentric offset field of width 10' and starting at r = rii,„ + 2'. The LF of the offset field was subtracted, bin by 
bin, from the cluster LF, taking the area proportion into account. The resulting LF was converted into an MF 
using the respective isochrone. It is woth noting that the cluster area, especially the core (r < Tc), may suffer 
significant incompletness due to stellar crowding. Therefore, our further results must be treated as a lower limit. 

The mass function 0(m) was approximated by a standard relation of the form log0(m) = — (1 + x) logm + foo 
where m is the stellar mass, x is the mass-function slope, and is a constant. The derived parameters allowed us 
to estimate the total mass Mtot and the total number of stars iVtot- These quantities were calculated extrapolating 
the MF from the main-sequence turn-off point down to the H-burning mass limit of 0.08 M© (see Maciejewski 
& Niedzielski 2007 and references therein for details). The mass function is plotted in Fig. [HI We obtained 
X = 1.7 ± 1.6, a value comparable within error bars to the universal initial mass function (IMF) ximf = 1.3 ±0.3 
given by Kroupa (2001). The high uncertainty is caused by a small range of covered stellar masses. The cluster 
was found to be in fact very rich and massive with A'tot = 31000 and Aftot = 12000 Mq. 

To describe a state of cluster dynamic evolution, the dynamical-evolution parameter r was calculated in the 
form r = where irciax is the relaxation time (see Maciejewski & Niedzielski 2007 for details). We derived 
logr = —0.3 which suggests that the cluster is dynamically younger than its relaxation time. 

4 Summary 

Berkeley 53 was found to be a rich and massive open cluster belonging to the Perseus Arm. Its age exceeds 
1 Gy, but it seems to be very young in the context of its dynamical evolution. A subtle discrepancy in the 
interstellar extinction law toward the cluster in the H and Ka, bands was detected. That suggests the presence 
of an abnormal interstellar extinction law toward the cluster. 
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^Following the Level 1 Requirement, according to Explanatory Supplement to the 2MASS All Sky Data Release and Extended 
Mission Products l|http://www.ipac.caltech.edu/2mass/releases/allsky/doc^ 





8 



ToK project MTKD-CT-2006-042514. This publication makes use of data products from the Two Micron All 

Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis 
Center /California Institute of Technology, funded by the National Aeronautics and Space Administration and 
the National Science Foundation. 

References 

[I] Dias, W.S., Alessi, B.S., Moitinho, A., Lepine, J.R.D.: 2002, A&A 389, 871 
[2] Giraldi, L., BerteUi, C, Bressan, A., et al.: 2002, A&A 391, 195 

[3] King, I.: 1966, AJ 71, 64 

[4] Kroupa, P.: 2001, MNRAS 322, 231 

[5] Maciejewski, G.: 2008, Acta Astron. 58, 389 

[6] Maciejewski, G., Niedzielski, A.: 2007, A&A 467, 1065 

[7] Majewski, S.R., Kron, R.G., Koo, D.C., Bershady, M.A.: 1994, PASP 106, 1258 
[8] Markov, H.: 2008, Bulgarian Astronomical Journal 10, 85 

[9] Mermilliod, J.C.: 1996, in The Origins, Evolution and Destinies of Binary Stars in Cluster, ed. E. F. Milone 
& J.-C. Mermilliod, ASP Conf. Ser. 90, 475 

[10] Niedzielski, A., Maciejewski, G., Czart, K.: 2003, Acta Astron. 53, 281 

[II] Odewahn, S.C., Bryja, C, Humphreys, R.M.: 1992, PASP 104, 553 

[12] Ruprecht, J.: 1966, BAICz 17, 33 

[13] Schlegel, D.J., Finkbeiner, D.P., Davis, M.: 1998, AJ 500, 525 

[14] Setteducati, A.F., Weaver, H.F.: 1960, in Newly Found Star Clusters, Radio Astronomy Laboratory, Berke- 
ley 

[15] Strutskie, M.F., Cutri, M.F., Stiening, R., et al.: 2006, AJ 131, 1163 



9 



